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Abstract 

An exact asymptotic analysis is performed for the two-point correlation func- 
tion C(r,t) in dissipative Burgers turbulence with bounded initial data, in 
arbitrary spatial dimension al. Contrary to the usual scaling hypothesis of a 
single dynamic length scale, it is found that C contains two dynamic scales: a 
diffusive scale Id ~ t 1//2 for very large r, and a super-diffusive scale L(t) ~ t a 
for r <C Id, where a = {d + l)/(d + 2). The consequences for conventional 
scaling theory are discussed. Finally, some simple scaling arguments are pre- 
sented within the 'toy model' of disordered systems theory, which may be 
exactly mapped onto the current problem. 
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I. INTRODUCTION 



Burgers equation (BE) has found many applications, in both theoretical and practical 
senses, over the years since its birth in 1940. It was originally proposed [Q to describe 
wave propagation in weakly dissipative media, and, in fact, it is now appreciated that 
within this large class of phenomena there are only two model descriptions in the limit of 
weakly non-linear waves; namely BE, and the Korteweg-de Vries equation [[J . In later years 
BE was scrutinised by the turbulence community as a simplified model of Navier-Stokes 
turbulence, and thus 'Burgulence' was conceived. The applications of BE were boosted 
again in 1986 when Kardar, Parisi and Zhang proposed that the BE with a stochastic 
source described the non-equilibrium evolution of a class of interface models [Q. Under a 
non-linear transformation, this noisy BE was seen to describe another rich class of systems, 
namely directed polymers in random media which have applications in wetting || , disordered 
magnets ||, and the pinning of flux lines in superconductors J7j. BE has also received 
attention as an approximate model for the formation of large scale structures in the universe 

II- 

Naturally, with such a wide range of physical applications, BE has attracted a great deal 
of theoretical attention. In the years subsequent to the revolution in critical phenomena - 
when the ideas of scaling and universality have become so prevalent ||TU[ - most theoretical 



ideas concerning BE are formulated within a 'scaling picture'. Although this is a convenient 
language for many phenomena, it must be realised that without a formal renormalization 
group (RG) description, scaling must be supported by strong physical insight and not merely 
'hand- waving' arguments. As an example, the physics of domain growth in quenched fer- 
romagnets has been very well understood on the basis of scaling arguments [|ll|], although 
no explicit RG calculations have been performed away from the critical temperature. The 
domain morphology of this problem provides an excellent basis for scaling since it is clear 
that the growing domain scale acts as a well-defined measure of dynamic correlations (with 
the caveat that scalar order parameter domain growth has more subtle scaling due to the 
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existence of a microscopic scale - the domain wall thickness). The concept of dynamical 
scaling is also supported in this field by a number of exact calculations (e.g. the one dimen- 
sional Glauber model ||12|| , and the large- n limit of the time- dependent Ginzburg-Landau 



equations [0]), and by a large number of numerical simulations. 

The existence of scaling is not so well established in BE, although most workers would 
agree that it is a convenient hypothesis, given the complexity of the problem. The analytic 
approach used by Burgers jl4|, and later by Kida |fL5| , certainly demonstrated the existence of 



an important length scale, which may be considered as the mean shock wave separation. The 
mute point is whether this is the single dominant length in the problem. If so, then one has 
dynamical scaling in its simplest form, and many quantities may be subsequently obtained 
by scaling arguments. What is lacking in the previous work on BE is an explicit solvable 
case in which scaling is seen to emerge cleanly. In order to achieve this it is necessary to 
calculate the form of some correlation function, which entails more difficulties than studying 
for instance the mean energy decay. Our intention here is to present such a calculation 
for a class of initial conditions in which the velocity potential is a bounded, discontinuous, 
random function. In this case exact calculations are possible and we may extract the form 
of the velocity-velocity correlation function E(r,t) for arbitrary spatial dimension d. We 
find that there exists dynamical scaling, but that it is controlled by two length scales, 
rather than one - a diffusive scale Id for large distances, and a super- diffusive scale L{t) for 
small distances. The details underlying this remark will be given below, but the important 
conceptual point is that if two length scales are playing a scaling role, then their ratio (which 
is, of course, dimensionless) may play a hidden role in subsequent scaling arguments. Thus, 
simple dimensional analysis is likely to fail. We shall see an explicit demonstration of this 
as we proceed. 

The outline of the paper is as follows. In the next section we introduce BE, discuss 
various choices of initial conditions, and briefly describe a few analytic steps which are 
required before the calculation proper. By adopting the initial condition mentioned above, 
along with some new analytic methods, we are able to calculate E(r, t), and we give explicit 



forms for its asymptotic behaviour for small, and large distances. This is presented in the 
lengthy section IV, section III being a warm-up exercise to calculate the mean energy decay 
(which has previously appeared in print ]16[). In section V we make a (formally exact) 
connection between BE and a popular 'toy model' in disordered systems theory. We then 
present simple scaling arguments for the toy model which yield partial agreement with the 
more complicated scaling picture that emerges from sections III and IV. The paper ends 
with our conclusions. 



II. DEFINITION OF THE MODEL 

BE is a partial differential equation written in terms of a velocity field v(x, t): 

3 t v = zA7 2 v - (v ■ V)v . (1) 

The field is taken to be irrotational, which allows one to express the equation solely in terms 
of the velocity potential defined via v = — V0. Explicitly one has 

<9^ = z/V 2 0+(l/2)(V0) 2 . (2) 

The equation is most commonly discussed in one spatial dimension, in the spirit of its 
application to non-linear waves ||. However, the <i-dimensional generalization given above 
is the canonical choice. One is interested in the evolution of the velocity field from some 
given initial condition, in the limit of vanishing viscosity, i.e. u — ► 0. This leads to strongly 
non- linear behaviour, otherwise known as the strong turbulence limit. We shall make this 
limit more explicit in terms of a dynamic Reynolds number as we proceed. 

The initial conditions we shall study are random functions vo(x), and as such are defined 
in terms of a distribution function -P[vo]. Naturally, there is an enormous choice available 
for P. Burgers WM studied perhaps the most natural, namely a Gaussian distribution of 



velocities with 5-function correlations: 



P B [v ] ~ exp (-(1/2D) J d d x v 2 ) . (3) 



His analysis was confined to d = 1 where a controlled analytic study was possible. The main 
result to emerge was that the velocity field forms into shock waves separated by smooth 
regions, and that the shocks become more dilute as time proceeds; the mean shock wave 
separation increasing as L s ~ t 2//3 . Dimensional arguments indicate that above d — 2, 
the asymptotic properties are dominated by diffusion (i.e. the shock waves disappear and 
the field diffusively vanishes), so that the dominant length scale is then a diffusion scale 
growing as t 1//2 . In precisely two dimensions |17], diffusion is still the dominant process, but 



logarithmic corrections are expected for quantities such as the mean energy decay S(t) = (v 2 ) 
(where here and henceforth, angled brackets indicate an average over the ensemble of initial 
conditions). Generalizations of the Gaussian form of the initial conditions (for example, 
defining different power spectra in Fourier space) have been studied previously |8|, |15| , |l8| , [19| . 
and scaling arguments have provided a broad classification for the time-dependence of the 
length scale L s (t) 

One may also consider initial distributions in terms of the velocity potential. A partic- 
ular class of these is to divide the system into cells of size l d , and to assign a value of 0o 
independently within each cell [|T^,|T^]. In this case we may write the distribution as 

P[V ] = II Pi^cell) ■ (4) 
cells 

It is important to distinguish between cell distributions p which allow bounded or unbounded 
values of their argument. We shall see that distributions of the former class (such as a 
top-hat function, or a cut-off exponential distribution) constitute a particular universality 
class, whereas those of the latter (such as a Gaussian or a power-law distribution) have 
different scaling properties. In the present work we shall be interested solely in random 
initial conditions of the former type, by demanding the cell distribution function to be 
defined only for a finite range of the velocity potential. Furthermore, one may show that all 
such distributions lead to the same asymptotic behaviour when the width of the distribution 
is large (but still finite), and we therefore concentrate on the simplest case, namely a top- 
hat function. [This is strictly true only for distributions which fall to zero discontinuously. 



Distributions which vanish at some finite value of their argument in a smooth manner require 
a separate analysis.] Explicitly we choose 

fl($-|0ol) 

where 8(z) is the Heaviside unit function 



Analytic progress has been possible in BE over the years, since for a given initial con- 



dition, one may exactly integrate the equation. This is due to the Hopf-Cole pl|j2^] trans- 
formation which linearizes BE. By defining u>(x, t) = exp[<p(x.,t)/2v] and substituting into 
eq.(|2]), one may see that w satisfies the linear diffusion equation which is immediately solved 
in terms of the heat kernel g{x,t) = (Aiiut)~ d ^ 2 exp[— x 2 /Avt]. Re-expressing the solution in 
terms of the velocity potential, one has the explicit solution of BE in the form 

0(x,t) = 2z/ln||A^(x-y,t)exp[0o(y)/2^]} . (6) 

The main analytic effort is to now perform averages over the initial distribution P. 
We shall accomplish this by making the following integral representation of the logarithm 
function in the above expression: 



oo 



lu(-) = / ^( e ~ u -e- uz ) . (7) 



o 



This representation has proven useful [23| in calculations in disordered systems theory as 



an alternative to the replica method, and has also been used previously in problems related 



to BE [|l(],|T7],^4]] ) . We therefore have the solution of BE in the form 

du , 



{x,t) = 2vj^[e- u -^(u,x,t)}, (8) 



where 



ip = expi^-ujd d y #(x- y, t) exp [0 o (y)/2z/]| . (9) 

Our main focus in this work is to calculate the velocity- velocity correlation function 
defined by 



£(r,t) = (l/2)(v(r,t)-v((M)> (10) 

which may be easily related (with the aid of translational invariance) to the two-point 
correlation function for the velocity potential via 

£(r,t) = (l/4)V 2 C(r,t) , (11) 

where 

C(r,t) = ([0(r,t)-0(O,t)] 2 ) . (12) 

The mean energy decay is given simply by £(t) = E(0,t), and we shall present a con- 
densed derivation of this quantity in the next section, before tackling the much harder task 
of calculating C(r, £). Results for E(0,t) have been presented before fL6f , but it is useful to 



sketch the derivation here in order to set up the necessary formalism required in section IV, 
along with revealing the important time scale in the problem. 

III. CALCULATION OF THE ENERGY DECAY 



In previous studies of BE p jl^jlql , it is more common to infer the energy decay from 
a scaling argument, once one has calculated the important dynamic length scale L s (t). On 
dimensional grounds one would like to infer £ ~ L 2 s /t 2 . This scaling relation certainly holds 
true in many situations, but it is by no means a universal result. We shall attack the problem 
from the opposite direction, by first calculating the energy decay explicitly. We shall then 
read off important length scales from the correlation function in section IV - comparing 
the two independent results will then allow us to see if dimensional analysis holds for our 
particular choice of the initial distribution. 

Expressing £ in terms of the velocity potential, one may see from averaging eq.(0) over 
the initial distribution, that 

£ = d t {<f>(x,t)) . (13) 
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So in order to determine the energy decay, we need only calculate the mean velocity potential, 
which in turn is related to the average of the function ip from eq.(Q). In fact a very similar 
function will be central in the evaluation of the correlation function, so it is useful to dedicate 
a few lines to deriving an explicit expression for (ijj) . 

In order to perform the average it is necessary to impose a lattice scale a - this is because 
generally the initial condition average has the form of a functional integral, which is only 
strictly defined on a lattice. We shall find that for all but the shortest times (set by to - 
the time for diffusion over the cell size 1) this scale a disappears from all physical quantities, 
and is replaced by the cell scale I which defines the correlation scale of the initial conditions. 
Explicitly we define a diffusion length = (4z/t) 1 / 2 and work in the limit a < I <^ lr>. 
In other words, the spatial smearing of the heat kernel is much greater than the cell size. 
Performing the initial condition average over ip using the distribution defined by eqs. (|4]) and 
(§), we find 

ln(V) =T d J Aln{-^- [Ei (ul d g(y,t)e- K °/ 2 ) - E x (vl d g(y,t)e*>' 2 )]} , (14) 

where E\(z) is the exponential integral p5[ and we have defined K to be the effective 
Reynolds number at zero time: i.e. Kq = (typical velocity x typical length)/// = $/za 
(We refer the reader to Appendix A where the initial condition average is performed ex- 
plicitly.) We may simplify this expression in two steps. Firstly we make the rescaling 
u — ► u{Tx l / 2 lr)/l) d e~ Ko / 2 and change the integration variable to s — y 2 /Aut. Secondly, we 
impose the strong turbulence limit by taking K ^> 1. This leads us to 



{^ 2 i D /iy 

K T(d/2 + 1) 



ln(V> = - Jtv^/o I \, L */M + 0(1/ 'K 2 ) , (15) 



where T(z) is the gamma function p5[ , and 



DC 

L p (u) = J ds s p [l - exp(-we~ s )] (16) 
o 

We refer the reader to Appendix B where it is shown that the integral may be evaluated for 
both small and large values of u with the result 
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uT(p+ 1) [1 - 2-( p+2 )w + 0(m 2 )] , u<l 
L p (u) = <( (17) 

lifgf + 7 pn(u)]* + 0([ln(«)]^ 1 ) , u » 1 . 

In order to find the mean velocity potential we must perform the w-integral as given 
in eq.(||). One may see that the w-integral is dominated by m < 1 (^> 1) when the ratio 
(Id/I) / K ^> 1 (<C 1) The former case occurs for very large times, and on performing the 
li-integral one obtains a diffusion result - i.e. £ ~ t~( d / 2+1 ). So, there exists a cross-over 
time t c ~ (I 2 /v)K^ d beyond which the non-linearity is irrelevant and the velocity potential 
relaxes according to diffusion. By taking the initial Reynolds number to be arbitrarily 
large, we may push t c to arbitrarily late times. The interesting non-linear regime occurs for 
to t <C t c in which case one must perform the u-integral using the large-w asymptotic 
form for (ip). In this case one finds (using the variable change a = [ln(w)] d//2+1 , and imposing 
a lower cut-off of 0(1) to the ^-integral) 

e = c t Wf™ ~ t ~, (is, 

where a = 2(d+ l)/(d+ 2) and Cd is a complicated ^-dependent constant. 

We may reinterpret this expression by defining a time-dependent Reynolds number K(t). 
For a typical velocity we take the square root of the mean energy decay, and for a typical 
(large) length scale we take Id (which we will justify a posteriori in section IV). Then we 
have 



K{t) 



C' d K 

\I'D 



l/(d+2) 



(19) 



where the constant C' d = n~ d / 2 T(d/2 + 2) has been chosen for future convenience. We can 
see that Kit) decays from its initial (very) large value with the power law t- d / 2 ( d + 2 ) until 
it becomes of order unity when t ~ t c . In this non-linear regime (defined by K(t) ^> 1) we 
may write the energy decay in the form 

£ ~ | K{t) 2 , (20) 
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which is cast into the form 'expected' from dimensional analysis, except that the dimension- 
less (but time dependent) Reynolds number is also present, which invalidates the prediction 
for the time dependence of S from dimensional considerations alone. 

The introduction of the time-dependent Reynolds number is useful, but must be justified 
by independently proving that lp is the typical (large) length scale in the non-linear regime. 
Alternatively one could insist on the dimensional prediction, in which case one would infer 
the important length scale to be L s ~ loK(t). (In fact, we shall see that the dynamic length 
scale ~ Id I ' K{t).) To place these results and speculations into a proper context one is forced 
to evaluate the scaling properties of the correlation function, which is a much more difficult 
task than the calculation of the mean energy decay. 



This section constitutes the heart of the paper in that we present the exact asymptotic 
forms for the correlation function C(r,t) in the non-linear regime to t <C t c . Unfortu- 
nately, in order to arrive at the required result, one must wade through a very long and 
technical calculation. So as not to burden the reader with details, all technical remarks will 
be relegated to appendices, with only the general flow of the analysis described in the main 
text. Firstly, we shall derive a general expression for C in the non- linear regime. In the 
subsequent sub-sections, we shall then analyse the asymptotic properties of C in the limits of 
r ^ Id and r ^> Id- As hinted at before, the main result of this analysis is the emergence of 
a new length scale which describes the small distance behaviour of the correlation function. 

It is convenient to define C in a symmetric way (cf. eq. (|T2"D ) : 



C(r, t) = M-v/2, t) - 0(r/2, t)] 2 ) = 2(0(0, tf) - 2(0(-r/2, t)0(r/2, t)) . (21) 



By utilizing the integral representation of the logarithm function twice, we may rewrite the 
bilinear combinations of velocity potentials in terms of integrals. This yields 



IV. CALCULATION OF THE CORRELATION FUNCTION 




(22) 



o 



o 
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where 



^(u, v, r, t) = (exp 



J d d y {ug{y - r/2, t) + vg(y + r/2, t)) e 



0o(y)/2^ 



(23) 



and we have utilized the property of translational invariance. 

In an analogous fashion to the averaging performed in the previous section, the average 
over the initial conditions may be performed in a straight-forward manner (see Appendix 
A), yielding a complicated expression for \l/ in terms of the exponential integral function. 
However, great simplification may be made by taking the limit Kq ^> 1. In this case we 
reduce \I/ to the form 



ln[$(u, v, r, £)] = —el(u, v, R) + 0{l/K { 



(24) 



where e = T(d/2 + 2)/K(t) d+2 < 1, and 
I(u, v, R 



2 r d ( 2 y-R 

rf 2/ 12/ ^~ 



we 



-(y-R/2) 2 



f e 



-(y+R/2) 2 



-(y-R/2)2 + t , e -(y+R/2)2 



we 



x < 1 — exp —we 



,-(y-R/2) 2 _ , ue -(y+R/2) 3 



(25) 



and we have defined R = r//^. 

At this point of the discussion it is convenient to consider the small and large distance 
behaviours of C separately. 



A. Small distance scaling 

To ascertain the small distance properties of C we need to perturbatively evaluate the 
above integrals in a power series in R 1. Although one may attempt this directly on the 
form of the integrals as given by eqs.(p3|) and (p5|), it is far more efficient to transform the 
function / beforehand into a natural power series in R 2 . The procedure for this is described 
in Appendix C, with the result 

oo 

I(u,v,R) = Y,(R 2 Y F p (u,v) , (26) 

p=0 
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where the functions F p have the integral form 

( — uyQ d )^ 

v) = r( P + i)r(p + rf/2 + i) L ^ u + v) ■ (27) 

We are interested in the non-linear regime K(t) ^> 1, and in this case the (u,v) integrals 
are dominated by u ^> 1 and v 1. Therefore we expand the integral L p+ d/2 appearing in 
eq. fl27D in powers of A = \n(u + v) ^> 1 (see Appendix B). The function F p may now be 
expressed as 

Fp{U,V) = T(p + l)T(p + d/2 + l) 



(28) 

where we have defined 

X P (u,v;q) = (uvd u d v ) p A q 

= f p (u, v- q)A q ~ 1 + g p (u, v- q)A q ~ 2 + 0(A^ 3 ) . (29) 

More details of these steps, along with the explicit form of the coefficients {f p } and {g p } 
may be found in Appendix D. 

Now that we have a workable series for I in the non-linear regime, it is possible to 
expand ^(u, v,r, t) in powers of R 2 , such that the coefficients are various combinations of 
the functions F p . The integrals over u and v may then be performed (see Appendix E for 
details), and one has the final result 

8z/ 2 (d + 2) \d + 2j w 



3{d + 2) 2 \d + 2J w 45(d + 2) 2 (d + 4) U + 2/ W 



4{d + 5){d+7) r fd + 8^ 



K{tf R 8 + 0(R m 



63(d + 2) 2 (d + 4) 2 (d + 6) V rf + 2 , 
+••• (30) 
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Several remarks are now in order. Firstly, the above result is given (after much effort) 
to quite high order in R 2 . One is obliged to do this to unambiguously determine the scaling 
properties of the correlation function. Secondly, the result for C has been written in such a 
way as to stress the form of the scaling. It turns out that the dominant term at each order 
of R 2 vanishes exactly, except for the dominant term at order R 2 - this explains why this 
term stands alone in the above expression. The sub-dominant terms from each order are 
non-zero, and are grouped together within the braces. The remaining terms play no role in 
determining the leading scaling behaviour and are indicated by the ellipsis. The fact that 
the dominant terms vanish means that the leading R 2 term can play no part in the scaling 
form of the correlation function. However, the terms in braces have a natural scaling form 
which allows us to read off a dynamic length scale. Explicitly we may recast the above 
expression into the scaling form (neglecting constants) 

L(fi) + fe) S {W), 
where S(x) is the scaling function, and the dynamic length scale is L(t) = ljj/K(t) ~ t a 
with a = (d + l)/(d + 2). 

We see that in the non-linear regime, the dynamic length scale is much smaller than the 
diffusive scale Id, although it is growing faster. This gives us another view of the cross-over 
from non-linear to linear evolution; i.e. the dynamic Reynolds number becomes of order 
unity when the scale L(t) becomes of the same order as Id- 

As a final remark in this section, we may obtain the velocity-velocity correlation function 
from eq. ([31]) with the use of eq. (|TT|) . One finds 



C(r,t) ~ v 2 



(31) 



v \ 2 x I r 



E { r,t) ~ Sit) + {-) S{ W) ) . (32) 

Again, it is interesting to see that the mean energy decay £ (t) is not part of the scaling 
function, which explains the difficulties encountered in the previous section with simple 
dimensional analysis. It remains to show the scaling importance of the diffusive scale - this 
will be accomplished in the next section. 
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B. Large distance scaling 



The scaling form for the correlation function for very large |r| may be obtained with 
relatively little effort. Starting with C(r,t) expressed in terms of the function I(u, f,R) as 
given in eqs.([22|), ([24]) and ([25]), we may express I by the following series (cf. Appendix B): 

00 



/(«,«, R) = J2 n ' W2+1) E C n m u m v n ~ m exp 



71 ' 

n=l m=0 



m r 

— (n — m)R 
n 



(33) 



As C(r, t) is non-zero for |r| — ► 00, it is convenient to measure correlations with respect to 
the asymptotic value C(oo,£). Thus we define 



<5C = CYoo, t) - C(r, t) = 8u 2 f— I —mu, v, r, t) - u, 00, t)] 

J U J V 


= 8z/ 2 y -^{exp[-eI(u,v,R)] - exp[-e/(u, -u, 00)]} . (34) 







From eq.fl3"3]) it is easy to see 



I(u , m) = g „- (u . + „.) = , (35) 

where we have rewritten the sum in terms of the familiar integral (using the integral 
representation shown in eq. ( |C6|) ). In the non-linear regime, we are interested in the large 
(u, v) behaviour, which according to eq. (17) gives us 

I(u,v,oa) ~ - 2) {[Hu)f 2+1 + {Mv)f^} . (36) 

Returning to eq.(|33|) the large- |R| form for / may be written as 



I(u, v, R) = I(u, v, 00) - 2~W 2+1 ) uve- R2/2 + ■■■ . (37) 



Substituting this result into eq.fl34) gives the leading term of SC as 



SC ~ ^ J (K? e~ R2 ' 2 , (38) 



where 
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oo 




(39) 



c 



(c being a number of order unity). 

The integral may be performed by steepest descents in the non-linear regime (details 
in Appendix B) with the result that 5C has the final asymptotic form (neglecting overall 
constants) 



where C'^ = d(d/2 + l)~( 1+2 / d ). We see from this expression that the diffusion scale Id is 
the natural scaling length for the correlation function, for very large distances. One may 
ascertain the range of validity of the above expression by calculating the contribution from 



Before ending this rather technical section, we shall recap the main results obtained. By 
performing an exact analysis on the correlation function C(r, t) in the non-linear regime, we 
have been able to confirm that there indeed exists dynamical scaling, albeit of a rather subtle 
type. The small distance properties of C are governed by a scaling length L{t) ~ lo/K(t), 
but the dominant term in C is singular - i.e. it may not be included into the scaling form. 
This indicates why the form of the mean energy decay found in the previous section was not 
obtainable by simple dimensional analysis. The scale L(t) is much smaller than the diffusion 
scale, but grows faster - the non-linear regime crosses over to simple diffusion when these 
two length scales become compatible. The large distance scaling was found to be more 
conventional, in that the dominant part of C (with respect to its asymptotic value) is a 
simple function of r/l D , albeit with a complicated prefactor; thus indicating that l D acts as 
the dynamic scale for the correlation function at very large distances. 



5C ~ u 2 K(tY A ~ d2 



exp exp[— r /2l D ] 



(40) 
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V. SCALING ARGUMENTS 



This section has two purposes. Firstly, we shall show an exact mapping between the 



solution of BE, and the free energy of a popular toy model |26| in the field of disordered 
systems. Secondly, we shall perform some simple scaling calculations p7j on the latter model 
to extract the form of the mean energy decay S(t) in the original BE problem. These scaling 



calculations are very similar in spirit to the original calculations of Burgers []14j and Kida 



The toy model in question is simply described. Consider a single particle in a potential 
composed of a harmonic background plus a random potential V(x). If the particle is in 
contact with a thermal reservoir, we may write the partition function for the particle as 



J d d x exp ^x 2 + V(x) | , (41) 



\2ti J 

where (3 is the inverse temperature, and we have normalized Z with respect to the harmonic 
background. For a given realization of the disorder potential V, we may calculate the free 
energy of the particle from Fy = —{1/0) ln(Z). 

At this stage we compare the expression for the free energy, with the exact solution of BE 
(evaluated at the origin) as given by eq.(|6]). We see that there exists an exact equivalence 
between the two quantities, if one makes the following connections: 0(0, t) «-> —Fy, 2v <-> 
1/(3, t <r+ l//i, and 0o( x ) ^ — ^( x )- This correspondence holds regardless of the particular 
distribution of initial conditions (or equivalently, disorder). 

To proceed with the description of the toy model, two quantities one may be interested 
in calculating are the quenched free energy F = (Fy), and the mean square displacement of 
the particle 

( x 2 ) = (^J d d x x 2 exp {-/? \^x 2 + V(x)] Yf 

= (-(2/(3)8, \n(Z)} = 2d,F . (42) 



By utilizing the correspondence with BE, we may relate the mean square displacement to 
the quantity S(t) in BE. Explicitly we write 
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(x 2 ) = 2d^F = -2d 1/t (<f>(0, £)} 

= 2t 2 d t ((j ) (0,t)) = 2t 2 £(t) , (43) 

where we have made use of eq.([U|). So we have been able to show that the dimensional 
prediction for the mean energy decay, namely £(t) ~ L s (t) 2 /t 2 , has a formal interpretation 
in terms of the toy model, so long as we interpret L s (t) as the root-mean-square displacement 
of the particle. 

We shall now derive an approximate form for L s (t) within the toy model formulation. 
Consider first, the top-hat distribution that has been the subject of the present work. We 
take P[V] = n p(V C eii), with p(V) = 0(D - \V\)/2D. The strong turbulence limit of BE 

cells 

corresponds to the zero temperature limit of the toy model. In this case, the particle will 
be trapped in the lowest potential energy minimum, within a given realization. In this case 
we may estimate the excursion of the particle by calculating the probability q(r, U*) for the 
lowest potential site to be located at a distance r from the origin, and to have an energy 
U*. This will be proportional to the probability that all sites within a radius r of the origin 
have a potential energy higher than U*. 

For a general potential distribution p(V), we may write 

oo 

q(r,U*)~p(U*-iir 2 /2) J[ J dV x p(V x ) . (44) 

l x l< r U*-iu?/2 

In the BE analogy we are interested in long times, so within the toy model we need to take 
H to be small, i.e., the harmonic background is taken to be very 'flat'. 

Restricting our attention to the top hat distribution, we see that a flat harmonic back- 
ground implies that the minimal energy U* will be close to the lower bound of the random 
potential —D. We therefore set U* = —D + SU, where \SU\ D. The above expression 
then reduces to 



rtr.m-jL n 

|x|<r 



SU fix' 



(45) 



2D AD 

It is now straightforward to exponentiate the term in square brackets, so as to transform the 
product over x as a spatial integral in the exponential. Taking fir 2 D (which is justified 
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a posteriori by the form of r typ given below) allows the integral to be performed simply, 
yielding the final result 



V( r , 5) ~ jjj ex P 



5 fr\ d , , s fir d+2 



(46) 



D \lj v ' l d D 

where C\ and C2 are constants. For this distribution of potential minima at distance r from 
the origin, we can read off a scaling relation between the typical value of r and /x; namely 

r w ~ (-) • (47) 

Making the correspondence with BE, we identify r typ <-> L s (i) and p,/D (I^Kq) . Com- 
bining the above result with eq . (|43|) , we see that we have derived the correct form of the 
mean energy decay as calculated previously in section III, cf. eq.( [i~8|) ; although the precise 
value of the prefactor may not be obtained by this simple scaling argument. 

The above result may also be cast into the form L s (t) ~ lr>K(t). This is guaranteed under 
the scaling hypothesis, but what is interesting is that such a length scale plays no role in the 
actual dynamical scaling as defined by the behaviour of the two-point correlation function. 
In other words, although we may calculate L s (or rather rtyp) from scaling considerations of 
the toy model, this length scale is not a dynamic scaling length - for instance, it could not 
be used to collapse the correlation function in a scaling plot. 

To end this section, we mention that the toy model may be analysed for other types 
of distribution. If one takes the disorder distribution to be of the cellular type, with p a 
gaussian, then one may rederive the result of Kida Hl5f ; namely that L s (t) is a diffusive 
scale up to logarithmic corrections. Alternatively one may consider the toy model in d = 1 
with a disorder distribution corresponding to the original Burgers choice, namely P[V] ~ 
exp[— J dx (dV/dx) 2 }. This particular scenario has been studied in detail recently, using a 
new replica approach ||28|| . Although the essential Burgers scaling (L s (t) ~ £ 2 / 3 ) is easily 
recovered, the calculation of prefactors is more difficult. It is an important test of various 
approaches as to whether they can quantitatively predict the correct prefactor. As far as we 
are aware, this is still an open problem, although there are a number of approximate values 
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VI. CONCLUSIONS 



This paper has been concerned with proving dynamical scaling for Burgers equation with 
random initial conditions. Exact calculations have been possible for a distribution of the 
initial velocity potential which has a large but finite region of support, and is uncorrelated 
in space. In section III we calculated the mean energy decay £ (t) in the non-linear regime 
(i.e. in the temporal regime in which the Reynolds number K(t) ^> 1). Explicitly we found 

£ ~ | K(t) 2 . (48) 

Dimensional arguments applied at this stage would then predict that there exists a dynamical 
length scale L s {t) ~ l^Kit), where Id ~ {yt) 1 ! 2 is the diffusion scale. 

In section IV we set out to establish this result by calculating the two-point correlation 
function C(r,t). This task is non-trivial, but we were able to extract the small and large 
scale asymptotics of C. There were two unexpected results. Firstly, the small and large scale 
forms of C, although assuming a scaling form, have different dynamic length scales. For the 
small distance scaling, the dynamic length scale was found to be L(t) ~ lu/K(t), and the 
dominant part of C in this spatial regime is singular, and can not be included in the scaling 
function. In terms of the velocity- velocity correlation function E(r,t), this singular piece 
is exactly £(t), with the scaling part of the correlation function describing the non-local 
properties of E; i.e. 

E< -^ ~ £ < () + (B ' s im) ' (49) 

The scaling function S(x) has a power series expansion in x 2 , the first four coefficients of 
which may be inferred from eq . (|30"D . The large distance scaling was found to be described 
by the diffusive scale Id which confirmed a posteriori the choice of this length scale in 
constructing the dynamic Reynolds number. The leading term in this large distance regime 
was found to be 

8C ~ v 2 K{t) {i ~ d2)/d eMCdK{t) 2{d+2)/d ] exp[-r 2 /2/^] , (50) 
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The second unexpected result is that neither of the two dynamic length scales coincide with 
the scale L s found from dimensional considerations of S(t). This is explained in part by 
the fact that the local energy decay is singular, and not contained in the scaling form for 
E. This result is similar to the cases in critical phenomena where caution is required in 
reading off scaling dimensions from composite operators (like v 2 ) [|Tl]]. Generally, these 
composite operators have their own scaling dimension, which may be related to the scaling 
of a two-point correlation function only via a small distance expansion (otherwise known as 
an operator product expansion.) 

In section V we introduced a mapping between BE and a toy model which is well known 
from the field of disordered systems - that being a thermal particle in a harmonic background 
along with a random potential. By considering simple scaling calculations on the toy model, 
we were able to reproduce the previous result for £ as found in section III (up to prefac- 
tors). The scale L s appears naturally in the toy model (as the typical root-mean-square 
displacement of the particle), and it is interesting that this scale is not a true dynamical 
scaling length, in that it may not be used to collapse the two-point correlation function in a 
scaling plot. It would be interesting to understand this result more fully, either by perform- 
ing similar calculations to those described here for other choices of initial distributions, or 
by calculating higher-order correlation functions, to see if more exotic forms of scaling (like 
multiscaling or intermittency) are present in these simple models. The connection between 
the toy model and BE may be of some mutual aid in both fields, at least in supplying new 
intuitive understanding of these complementary problems. 

The author would like to thank Sergei Esipov for important discussions in the early 
stages of this investigation, and Alan Bray for an enlightening conversation. The author ac- 
knowledges financial support from the Engineering and Physical Sciences Research Council. 
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APPENDIX A: 



In this appendix we outline the initial condition average over the function 

ij; = exp{-ujd*y <?(y,i)exp[0 o (y)/2i/]} , (Al) 

where we have taken advantage of translational invariance and set x = 0. (In fact, the 
translational invariance holds only for times greater than to, since clearly the cellular initial 
conditions allow only invariance under discrete transformations over a period I. Once the 
heat kernel has diffused beyond the cell scale, the continuous translational invariance is 
recovered.) The initial distribution is as described in section II, namely, we divide space into 
cells of volume l d , and within each cell we assign the velocity potential to be an independent 
random number drawn from a top hat distribution of width 2$. Explicitly we have 

(V) = / 2><Myi) P[0oM0o] 

X - J d<MY) exp |- W exp[0 o (Y)/2H £ ^(y,t)| , (A2) 



5 2$ 



where Y labels the cells, and y$ labels the discretized points (on a scale of the lattice cut-off 
a) within a given cell. The integrals are easily performed in terms of the exponential integral 
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— / rf0oexp[-Aexp[0o/2z/]l = — \E x (Atr K °l 2 ) - E^Ae^ 2 )} , (A3) 

— $ 

where K = $/z/ is the initial Reynolds number as defined in the text. To reach the result 
shown in the text, namely eq.(rj), two more steps are required. Firstly the above result is 



re-exponentiated, so that the product over cells in eq. ( |A2| ) may be written in the exponent 
as a sum over cells. Secondly, we use the fact that for times greater than to = I 2 / v, the heat 
kernel has smeared beyond the cell scale, so that the sum over cells (on a scale I) and points 
within cells (on a scale a) may be replaced once more by a continuum spatial integral. In 
this way one finally arrives at the result given in the main text. 
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APPENDIX B: 



This appendix is dedicated to the asymptotic evaluation of two integrals that appear in 
the main text. 



The first integral L p (u) is defined in eq. ( IB ) , and is used repeatedly in the present work. 
We have 

oo 

L p (u) = J ds s p [l - exp(-ue' s )} . (Bl) 
o 

We require both the small and large-u forms for this integral. The former is trivially obtained 
by expanding the integrand as a power series in u. The asymptotic expansion of this integral 
for u 3> 1 is less simple. For the precision required in the calculations, we need both the 
dominant and sub-dominant terms. To extract these we proceed as follows. We notice 
that for large-ti, the second factor in the integrand behaves very much like a step function 
centered at s = ln(u). Thus, the dominant term will arise from replacing this factor by 
6*(ln(u) — s), and the sub-dominant term will arise from finding the leading error made by 
this approximation. 

So, explicitly we write 

InO) oo 

L p (u) = J ds s p + J ds s p T(s) , (B2) 



where T(s) = [1 — exp(— ue~ s )] — 6(\n(u) — s). The first term gives the dominant contribution 
to the integral, which equals [ln(«)] p+1 /(f> +1). To extract the main contribution from the 
second term, we replace s p by [ln(-u)] p and perform the integral over T(s): 



ds s p T{s) = [ln(n)] p J ds T{s) + 0{[ln{u)] p - 1 ) 





oo ln(u) 

= [ln(u)] p < J ds [1 — exp(— ue~ s )] — J ds exp(— ue 

Jn(«) 



> + 0([\n(u)} p - 1 ) . (B3) 



The first integral in the braces may be evaluated simply (using the variable change q = e s ) 
to give l) n+1 /nn\ , whilst the same variable change reduces the second integral to 

22 



ln(u) 



J ds exp(-ue~ s ) = J ^ e~ q = E X {1) + 0{e~ u /u) . 



(B4) 



Using the series expansion for the exponential integral [25] , we may combine the two results 
to give 



ds s p T(s) = 7[ln(w)] p + 0(pn(u)] 



p— 1\ 



(B5) 



where 7 = 0.57722... is Euler's constant. This is the required result for the subdominant 
terms. 

The second integral appears in the evaluation of the large distance scaling in section IV. 
Referring to eq.(^) we need to evaluate 



00 

J{K) = J du exp{-K(t)-^[\n(u)} d/2+1 } 



(B6) 



in the non-linear regime, K(t) ^> 1. The constant c is a number of order unity, and arises 
since we have used the large-w form for the integrand, and so we must cut off the integral 
at the lower end. For notational convenience let us consider 



M b (N) = J du exp{-(l/iV)[ln(n)] b } , (B7) 

1 

for N ^> 1, with b > 1. Then we can retrieve the integral of interest from J(K) = 
M d/2+1 (K d+2 ). 

In order to cast the integral into a form suitable for steepest descents, we make the 
variable change x = jV -1 '( 6-1 ) ln(-u). We then have 



M b (N) = N l/{b - l) j dx exp 




(B8) 



This integral is easily performed by steepest descents to give (neglecting overall 6-dependent 
constants): 

"(6-1) /N\^ b -^ 



M b (N) ~ N 1 ' 2 ^ 



cxp 



b V b 

from which one may retrieve the form given in eq. 



(B9) 
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APPENDIX C: 



In this appendix we give details of the manipulation of I(u, v, R) into a series expansion 
for small and large R. As given in eq.(p5|) we have 



I(u, v, R) 



dll d / 2 



,d j 2 y R 

dy <y 



Me -(y-R/2) 2 _ we -(y+R/2) 2 -i 



ue 



-(y-R/2)2 + ue -(y+R/2)2 



x < 1 — exp 



(CI) 



As a first step, we expand the exponential term in the second factor of the integrand to 
obtain 



\n+l 



n=l 



-AJu,v,R) 



(C2) 



where 



Jd d y {y 

y R 



Me -(y-R/2) 2+t;e -(y + R/2) 2 



we 



-(y-R/2) 2 + we -(y+R/2)5 



n-1 r 



Me -(y-R/2) 2 _ ve -(y+R/2) 2 



(C3) 



This cumbersome expression may be simplified greatly by expanding the brackets as binomial 
series in powers of u and v, and then performing the gaussian integrals over y. One is then 
left with 



A(u, v, R) = n - W2+1) J2 C^u m v n ~ m exp 



m=0 



771 



(n — m)R 



n 



(C4) 



Combining eqs.flC2]) and (|C4[ ) reproduces the large- i? form given in eq.(^). 

In order to cast A into a small- R form, we expand the exponential terms in eq. (pjj) and 
then binomially resum the series in u and v. This leaves us with 



(uvd u d v ) p (u + v) 



(C5) 



p^o nP P ] 

We now substitute this expression back into eq.( p2|) and make the integral representation 
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n -(p+d/2+l) 



[r(p + d/2 + l)]" 1 J ds s p+d/2 e- ns . (C6) 



This allows the sum over n to be performed as that for a geometric series, and we are left 
with 

I(u, v, R) = r( p + !)r(p + d / 2 + 1) ( M ^^) PL p+rf/2^ + «) , (C7) 
as given by eqs.(p6|) and (|27| ) in the main text. 



APPENDIX D: 



This appendix will give details of the rewriting of the function F p (u,v) in moving from 



eq.fl2T|) to fl28|) in the main text. As given by eq.(^), we have defined 



7-1 / x (-uvd u d v ) p , . 

= T( P+ in P+ d/2 + i) w(«+«) • ( D1 ) 

Since we are interested in the non-linear regime, we may expand the integral L as shown in 
Appendix B. Namely, 

Ap+d/2+1 

L p+d/2 (u + v)= {p - — - + 7 A p+d / 2 + 0(A p+d/2_1 ) , (D2) 

where A = ln(u + v). Then we define the quantities Xp{ u i v i <?) = ( uv d u d v ) p A q which allow 
us to rewrite the above expression as 



F p (u,v) 



T(p + l)T(p + d/2 + l) 



(D3) 



which is of the form given in eq.(^) in the main text. 

The above steps are largely a matter of redefinition of various quantities. We must extract 
explicit forms for the functions Xp( u -> v 'i <?)• Clearly the zeroth function is just Xo( u , v , o) — 
A q . The subsequent functions may be written as 

X P (u,v;q) = {uvd u d v ) p A q 

= f p (u, v; q)A q - 1 + g p (u, v; q)A q ' 2 + 0(A q - 3 ) . (D4) 

25 



In order to unambiguously determine the scaling properties of the correlation function, it 
is sufficient to explicitly evaluate C to 0(R 6 ). However, given the singular nature of the 
scaling in this problem, we shall proceed to calculate the 0(r 8 ) terms as well, as a useful 
check. This in turn necessitates calculating the coefficients f p and g p for p = 1,2, 3, 4. With 
the application of brute force algebra, we obtain: 

fi{u,v;q) = - q 



HI) 

f 2 (u, v;q)= + q- — ■ — - {u 2 - Auv + v 2 ) 

[U + Vf 
WV 

fz(u, v;q)= — q- — (u 4 — 26u 3 v + 66u 2 v 2 — 26uv 3 + v 4 ) 

(u + v) b 

111) 

U(u, v;q)= + q- (u 6 - 120u 5 v + 1191u 4 w 2 - 2A16u 3 v 3 + 1191«V - 120uv 5 + v e ) 



(D5) 



and 

gx(u,v;q) = - (q - 1)^ 

u 2 v 2 

g 2 (u,v;q)= - 2(q - l)/ a - q(q - I)- — ■ — r- 

[U + V J 4 

u 2 v 2 

g 3 (u, v- q)= - 3(q - l)/ 3 + q(q - 1) (17u 2 - 52uv + 17v 2 ) 

(u + v)° 

g 4 (u, v;q)= - 4(g - l)/ 4 - q{q - ^ (129m 4 - 16A8u 3 v + 3538u 2 v 2 - 1648m; 3 + 129i; 4 ) . 

(D6) 



APPENDIX E: 

In this final appendix we give a brief description of the final steps required in order to 
perform the u and v integrals, so as to derive the final form of the correlation function given 
in eq.fljg). 

Referring to the main text, we see that combining eqs.(p2"D, (El) and (p6|), we may write 
the correlation function in the form 
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oo , oo 



C(r,t) = 8u 2 I — I — exp[-eF (u, v)\ \ 1 - exp 

J U J V 



(El) 



where c is a number of order unity, required simply to cut off the integrals at their lower 
limits, given we are using the asymptotic form for the integrand in the non-linear regime. 
The remaining steps are easy to describe although rather tedious to perform in practice. We 
expand the exponential in the last factor of the above integrand in powers of R, and then 
integrate over u and v. We shall explicitly demonstrate this for the dominant part of the 
0(R 2 ) term. Using the asymptotic forms of the functions Fq and Fx from Appendix D, we 
have 



r, OO OO 

sv e f du r dv 



° (R2) ~ W 2 + 2)r W 2 + 2) J 'J- «"< "> 2 + 6XP 



r(d/2 + 2) 

(E2) 



Remembering that e = T(d/2 + 2)/K(t) d+2 , and using the form of fi from the previous 
appendix, we may rewrite this term as 



oo oo 



0{R 2 ) term ~ 8u 2 K{ty^ J du J dv (u + v)- 2 A d ' 2+l exp [-K{t)-^ A d ^ 2+1 ] . (E3) 

c c 

Now the double integral has the form 

X OO OO OO OO 

du I dv (u + v)~ 2 A(u + v)= du dw w~ 2 A{w) ~ / — A(u) , (E4) 



u 



c c c u c 



where the last step was achieved using integration by parts on u. The final integral may be 
easily performed by substituting x = ln(u), thus yielding the first term in eq. (|30|) . 

The integrals required at higher orders in R 2 may all be performed by changing vari- 
ables from (u, v) to (u, w = u + v) and performing the appropriate number of integrations 
by parts on u. In general the 'boundary' terms do not vanish, but they are negligible in 
the limit of interest, namely K(t) ^> 1. (This is because the general form of the func- 
tion A(w), which appears as a spectator in these integral manipulations, is proportional 
to exp[— K(t)~( d+2 ^ [ln(u>)] d / 2+1 ], which decays faster than any power.) A final point worth 
mentioning is that the sub-dominant terms (i.e. those terms involving the coefficients g p ) are 
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required since the (u,v) integrals over the leading terms containing f p all vanish for p > 1. 
This leads to the unusual scaling form described in the text, in which the dominant term 
may not be cast as part of the scaling function; and explains our cautionary calculation of 
the 0(R 8 ) term, which indeed confirms this singular scaling, i.e. the integrals over f\ vanish 
exactly. 
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